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INTRODUCTION 

Early  in  the  design  phase,  helicopter  studies  require  a great  deal  of  aero- 
dynamic information,  such  as  surface  air  loads  and  flow  directions,  and  various  compo- 
nent interference  velocities.  This  information  is  most  advantageous  before  any  one 
configuration  Is  identified  as  a baseline  model  for  wind  tunnel  testing.  In  addition, 
aerodynamic  data  required  for  the  most  efficient  design,  such  as  various  rotoi  Tuselage- 
induced  interference  velocities,  are  often  costly  and  extremely  difficult  to  obtain  in  a 
wind  tunnel  test.  Theoretical  and  experimental  works  in  helicopter  aeiodynomics  ore 
presented  in  References  1 to  10. 

The  induced  interference  effects  are  of  particular  importance  since  the 
separated  flow  in  areas  of  high  interference  velocities  is  the  major  source  of  aero- 
dynamic drag.  The  interference  velocities  result  in  steep  adverse  pressure  gradients 
that  may  lead  to  further  separation.  Furthermore,  the  separated  wake  from  rotor/ 
fuselage  components  creates  additional  interference  due  to  wake  blockage,  Ttiese 
effects  are  seen  in  Figure  1.  The  interacting  flow  fields  of  the  pylon  ramp,  pylon, 
shaft,  hub,  and  blade  shanks  produce  a strong  diffuser  effect  that  leads  to  gross  separa- 
tion, To  minimize  the  separation  and,  therefore,  the  induced  drag,  it  is  necessary  to 
eliminate  the  peak  interference  velocities  by  gross  redistribution  of  cross-sectionol 
area  and  to  contour  the  localized  interference  regions  around  the  pylotr  to  avoid  the 
sudden  deceleration  of  flow  in  the  boundary  layer. 

In  the  design  phase,  it  Is  desired  to  be  able  to  investigate  analytically  the 
local  interference  regions  and  the  effects  of  geometry  variations  on  the  separatior> 
characteristics  of  the  design.  The  advent  of  high-speed  computers  has  allowed  the 
application  of  finite  element  aerodynamic  methods  to  complex  aerodynamic  conflgura- 
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OUASI-ONE  DIMENSIONAL  FLOW  "OIFFUSOR  EFFECT 

HIGH  LOCAL  STREAMLINE  ACCELERATKTN 


Figure  1.  The  Hub-Pylon  Interference  Problem. 

tions.  As  a result,  surface  airloads  may  be  obtained  quickly  and  efficiently.  Flow 
direction  and  magnitude  can  be  obtained,  aiding  in  the  identification  of  potential 
separated  flow  regions  to  minimize  the  effect  on  aircraft  drag. 

This  report  presents  the  results  of  a study  to  examine  the  total  effect  of  tt>e 
presence  of  the  rotor  in  its  various  operating  conditions  on  the  pressure  distribution 
and  drag  (induced)  of  the  helicopter  components  beneath  the  rotor,  A potential  flow 
aerodynamic  program  with  suitable  rotor  wake  representation  is  developed  to  predict 
the  separation  characteristics  of  arbitrary  three-dimensional  bodies  with  lifting  sur- 
faces in  yawed  flow.  In  particular,  the  effect  of  the  main  rotor  blade  wake,  blade 
shank  wake,  and  hub  woke  on  the  separation  of  the  flow  over  the  hub,  shaft,  pylon, 
and  fuselage  of  a lifting  helicopter  in  forward  flight  is  analyzed. 
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The  report  is  presented  in  two  volumes.  This  volume,  Volume  I,  Theoreti- 
cal Formulation,  gives  the  theoretical  background  for  the  computer  program.  Volume 
II,  Program  SHAPES  User's  Manual,  presents  a detailed  description  of  the  computer 
program . 

This  volume  begins  with  the  background  description  of  the  helicopter  rotor- 
fuselage  aerodynamic  interference.  The  analytical  approach  to  the  problem  for 
separated  flow  and  the  helicopter  rotor  wake  representation  are  also  discussed.  A 
theoretical  foriviulation  for  incompressible  aerodynamics  with  separation  is  given  in 
Section  3.  Fundamental  aspects  of  wake  dynamics  are  discussed  in  Section  The 
numerical  formulation  of  the  aerodynamic  method  for  use  in  the  computer  program  is 
given  in  Section  5.  Section  6 presents  the  computed  results.  Concluding  remarks 
on  the  present  study  and  recommendations  for  future  research  are  also  provided.  A 
list  of  references  supporting  the  various  discussions  is  also  included. 
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SECTION  2 

HELICOPTER  AERODYNAMICS 

TLo  liollcoptor  configuration  with  its  numarous rotating  parts  is  a complex 
aorodytramic  shape  with  intersecting  wakes  and  separated  flows.  Analysis  of  these 
configurations  requires  the  use  of  computer  progranrs  designed  for  use  orr  high-speed 
computers.  Programs  of  this  type  that  can  ho  applied  to  configuiations  with  aoiody- 
namic  interference  are  reviewed  in  Subsection  2.1  . One  ol  these,  Program  SOUSSA, 
is  discussed  in  more  detail  in  Subsection  2.2,  and  its  application  to  lielicopter  aero- 
dynamics is  introduced.  The  last  subsection  presents  a brief  discussion  of  the  repre- 
sentation of  the  helicopter  rotor  wake. 

2.1  AERODYNAMIC  INTERFERENCE 

The  classical  approach  for  evaluating  the  aerodynamic  pressure  on  interfer- 
ing lifting  surfaces  (wings  or  rotors)  is  the  lifting  surface  method  (Reference  11). 

Calculations  have  boon  normally  limited  to  wings  or  other  lifting  surfaces  separately, 
except  for  a limited  number  of  cases  in  which  the  vehicle  could  bo  considered  as  a 
slender  Ixsdy.  However,  interaction  of  the  flow  fields  of  lifting  surfaces  and  the 
fuselage  results  in  intorforonco  effects  (Reforonco  12)  that  cannot  lie  adequately 
accouirtod  for  wlion  examining  the  liftii>g  surface  and  the  body  separately.  Tlio  advent 
of  modern  high-speed  computers  has  prompted  the  recent  development  of  a series  of 
now  methods,  generally  called  finite  element  methods.  Those  methods  offer  both 
flexibility  and  efficiency  and  can  bo  applied  to  complex  configurations.  All  the 
methods  are  based  on  potential  aerodynamics  and  considoi  various  distributions  of 
sources,  doublets,  vortices,  pressure  panels,  etc.,  on  and/oi  inside  the  surface  of 
the  aircraft . 

j 

i; 
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An  eaily  work  on  the  flow  field  utout^d  three-dimensional  LK>dies  hy  Mess 
and  Smith  (Reference  13)  uses  constat\t  stieiryth  source -elements  to  solve  the  pushlem  ol 
steady  subsonic  flow  around  nonlifting  bodies.  This  metlu^d  has  been  extended  K>  lift-  I 

ing  bodies  ( References  14-16)  by  including  doublet  and  vortex  panels,  Woodwaid's  i 

method  for  steady  subsonic  and  supei sonic  flow  (Reteiences  17  and  18)  is  a difteieni  | 

approoch  that  uses  liftiirg  surface  elenrents  Kvi  the  lepiesentation  of  the  body.  Little 

work  has  been  doire  in  unsteady  flow,  for  oscillotory  subsonic  aeuKlynomics,  exten-  ■ 

s 

> 

sions  of  the  doublet-lattice  method  (Reference  19)  oie  obtoined  by  eithei  plocing  i 

unsteady  lifting  surface  elements  on  the  suiface  of  the  l>ody  or  by  using  the  metlwd  of 

images  combined  with  slender  body  theory  (Reference  20).  In  the  supersonic  range, 

conrplex  configuiations  are  analyzed  in  Refeiences  17,  21,  and  22.  The  pivsgiam 

noted  in  Referei^ce  7 is  a modification  of  the  pivvgram  by  Rubbei  t and  Saaiis  (see 

Reference  14),  which  ii>  turn  is  based  on  the  original  proyiom  by  Mess  and  Smith 

(Reference  13). 

Of  all  the  methods  meistioned  above,  luxse  is  sidtTciently  geneiol  h.'  oiur- 
lyze  unsteady  rotor -fuselage  Interference  in  compressible  flow  since  there  exists  no 
frame  of  reference  witls  respect  to  which  all  the  surfaces  (boiTy  and  rotiul  oie  fixed. 

) 

However,  a general  method  for  unsteady  con\pressible  aerosfyiHimics  arouiui  complex 

configurations  (where  the  surface  is  allowed  to  move  with  resfsect  to  the  home  of  | 

referetsce  os  necessary  for  the  present  atsalysisT  has  recetstly  been  develo|,'ed  (Refer-  | 

ences  23  and  24).  This  method  provides  o unified  theory  foi  onolyzing  steady  and  | 

i 

unsfeody,  subsonic  and  sifsersonic  aerodynomics  for  ciwifvlex  oircioft  configmathvns  j 

i 

and  is  the  only  method  that  allows  for  the  complex  unsteady  roK>i -Uiseloge  intei - 
actiofs  in  compressible  flow.  Various  develofsments  are  given  in  References  23  to  33. 

This  method  is  readily  opfslied  to  compressible  (subsoisic  and  sipei sonic)  I 

flow  and  is  the  only  known  method  for  analyzing,  for  insKince,  unsteosfy  supeisvvnic  j] 


^ ^ 

I 

t 

flow  around  non -zero-thickness  configurations  (References  27  and  29).  The  technique 

is  implemented  in  the  computer  program  called  SOUSSA  (Steady  Oscillatory  Unsteady  ! 

Subsonic  and  Supersonic  Aerodynamics).  In  this  program,  the  motion  of  the  surface 

with  respect  to  the  frame  of  reference  is  assumed  to  consist  of  small  osr.illotions  with 

arbitrary  variation  in  time.  The  surface  of  the  aircraft  and  its  woke  ore  divided  into 

small  quadrilateral  elements  that  are  approximated  by  a hyperboloidal  surface  defined 

by  four  corner  points.  In  this  process,  the  continuity  of  the  surface  is  maintained, 

although  discontinuities  In  the  slopes  are  introduced.  The  unknown  is  assumed  to  be 

constant  within  each  element  and,  therefore,  the  integral  equation  reduces  to  a system 

of  linear  algebraic  equations.  It  should  be  noted  that  the  general  formulation  includes 

fully  unsteady  flow.  Therefore,  the  problem  of  interference  of  the  fuselage  with  the 

rotor  is  a particular  case  of  the  general  theory  of  Reference  23. 

2.2  PROGRAM  SOUSSA  j 

Program  SOUSSA  is  based  on  radically  new  integral  equations  proposed  first 
by  Dr.  Luigi  Morino  (Reference  23)  for  lifting  non-rigid  bodies  having  arbitrary  shapes 
and  motions  (including,  for  instance,  the  combined  motion  of  blades-fuseloge  for  a 
helicopter  In  forward  flight  or  maneuvering,  with  elastic  deformation  Included).  The 

method  is  based  upon  the  classical  Green's  theorem  approach,  which  was  never  before  j 

used  in  computational  aerodynamics.  While  the  representatian  of  the  velocity  poten- 
tial is  still  in  terms  of  source  and  doublets  distribution,  the  specific  nature  of  the 
representation  (which  is  not  assumed,  but  obtained  from  the  Green's  function  theory*) 
has  many  advantages  that  make  the  method  superior  to  any  other  existing  one.  Only 
the  subsonic  case  is  described  here.  First,  the  integral  equation  does  not  imply  the 

*For  example,  for  steady  incompressible  flows,  one  obtains  a distribution  of  source  with 
known  Intensity  eaual  to  the  downwash  (prescribed  from  the  boundary  conditions)  and  a 
distribution  of  doublets  eaual  to  the  unknown  potential  itself.  This  is  obtained  by 
solving  the  integral  equation  (see  Section  3). 
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normal  derivoHve  of  fhe  integral  representation  (as  the  others  do)  and,  therefore, 
there  ore  no  singularities  in  the  integral  equation  (see,  for  example.  Reference  23). 

This  yields  programming  simplicity  with  the  scheme -convergence  being  very  fast. 
Therefore,  fewer  elements  are  required  and  thus  the  computer  times  are  very  low. 
Second,  the  vorticity  distributions  over  the  boundary  layer  and  the  wakes  are  included 
(as  equivalent  doublet  distributions)  in  the  limit  as  the  viscosity  goes  to  zero  (zero- 
thickness boundary  layer  and  wakes).  Third,  Program  SOUSSA  is  very  simple  to  use 
since  the  type  of  distribution  (source  doublet,  or  vortex  elements)  is  prescribed  by  the 
formulation  itself  and  not  left  to  the  experience  of  the  engineer  using  the  program. 
Finally,  the  methodology  is  applicable  to  the  general  compressible  unsteady  flow.  In 
particular,  the  effect  of  the  compressibility  and  unsteadiness  for  the  rotor-blade 
motion  can  be  included  in  an  exact  way.  No  other  existing  theory  can  accomplish 
this. 

The  classical  approach  is  weak  in  the  modeling  of  separated  flows  since,  in 
the  separated  wake  (i.e.,  the  region  between  the  dividing  line  and  the  body),  the 
flow  field  is  not  potential . The  presence  of  vorticity  is  apparent  for  the  point 
(Figure  2).  Therefore,  formulations  of  potential  flow  with  boundary  layer  (where 
the  vorticity  is  limited  to  a thin  layer  near  the  body)  are  not  adequate  to  describe  the 
complete  flow  field  since  vorticity  is  not  included  in  the  neighborhood  of  the  point  Q^. 

The  correct  formulation  for  the  description  of  flow  fields  in  the  presence  of 

vorticity  (rotational,  i.e.,  nonpotential  flows)  is  given  in  Reference  30.  A simplified 

version  of  this  general,  exact  formulation  is  used.  This  consists  in  assuming  that  the 

vorticity  distributed  in  the  separated  wake  can  be  replaced  with  the  vorticity  of  a 

zero -thickness  separation  wake*  with  the  addition  of  a single  isolated  vortex  located 

*The  separation  wake  is  defined  here  as  a zero -thickness  vortex  sheet  (i.e.,  doublet 
layer.  Reference  23)  emanating  from  the  separation  line  and  composed  of  the  dividing 
streamlines. 
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FLUID  ORIGINATING 


Figure  2.  Separated  Flow. 

at  Q^,  The  location  of  the  separation  line  and  the  geometry  of  the  zero -thickness 
wake  along  with  the  geometry  of  the  vortex  line  is  assumed,  based  upon  experimental 
results,  physical  intuition  and  engineering  experience. 

Program  SOUSSA  can  be  used  to  solve  the  problem  If  the  value  of  the 
strength  of  the  vortex  is  known.  This  value  may  be  obtained  by  iteration  by  imposing, 
at  a prescribed  point,  that  the  separation  line  is  a stagnation  line.  In  addition, 
various  geometries  for  the  wake  and  various  locations  for  the  separation  line  will  be 
analyzed  by  iteration  and  imposing  the  condition  that  the  flow  field  is  tangent  to  the 
wake  surface  and  that  the  separation  line  corresponds  to  a stagnation  line.  The 
agreement  with  experimental  results  will  be  used  as  final  criterion  for  assessing  the 
improvement  and/or  correctness  in  the  solution. 
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2.3 


ROTOR  WAKE  GEOMETRY 


f 


The  complex  woke  geometry  of  a helicopter  during  cruising  flight  is  shown 
in  Figure  3.  The  most  important  parameters  are  the  distribution  and  time  history  of 
blade  loading*,  which  determine  the  internal  structure  of  the  wake  as  well  as  the 
varying  size,  shape,  and  deflection  angle  of  the  wake,  which  in  turn  determine  the 
position  of  the  wake  elements  with  respect  to  components  of  the  aircraft  such  as  the 
fuselage. 


BLADE  LOADING 


Figure  3,  Helicopter  Rotor  Wake  Geometry. 


An  approximation  to  the  blade  loading  and  vortex  wake  of  the  helicopter  is 
made  as  follows.  The  wake  for  each  blade  is  assumed  to  be  a helicoidal  surface  of 
vorticity  whose  axis  is  at  a wake  deflection  angle  X to  the  free  stream  velocity,  con- 
tained within  a circular  cylinder  of  fixed  radius  whose  axis  is  at  an  angle  (i  + x)  with 
respect  to  the  rotor  tip-path  plane.  This  rotor  wake  model  is  shown  in  Figure  4. 

*|t  should  be  noted  that  the  blade  loading  is  automatically  obtained  in  Program 
SOUSSA  if  the  geometry  of  the  wake  is  assumed. 
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Figure  4.  Rotor  Woke  Model. 


Presently,  woke  controction  is  not  included.  Because  of  the  asymmetric  geometry  of 
the  rotor  wake  as  well  as  the  presence  of  wake  interactions  of  the  rotor  with  the 
separation  line  and  rotor/fuselage  wake  interaction,  simplifying  assumptions  are  made. 
These  are: 

e The  vortex  patterns  from  each  blade  do  not  intersect 

e The  rotor  woke  can  be  truncated  above  the  fuselage 

e The  rotor  wake  is  described  so  that  it  does  not  intersect 
the  separation  line  or  the  isolated  vorticity 

e The  separation  line  and  the  isolated  vorticity  are  not 
coplanar. 

It  is  shown  in  References  34  and  35  that  when  the  rotor  has  a smell  incidence 
angle  with  respect  to  the  freestream  (less  than  30°),  the  wake  rapidly  rolls  up  into  two 
concentrated  vortices  similar  to  the  tip  vortices  of  a monoplane  wing.  On  this  basis, 
the  wake  model  chosen  for  the  helicopter  is  equivalent  to  two  vortices  trailed  from  the 
ends  of  the  lateral  diameter  of  the  rotor  disk.  The  strength  of  the  vortices  is  obtained 
by  assuming  that  the  rotor  behaves  as  an  elliptical  monoplane  wing  of  span  2R  (the 
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Giauert  hypothesis)  whose  lift,  i.e.,  thrust,  T,  is  given  by 


■f [■  -(5)1' 


where 


= air  density 
= flight  velocity 

= midspan  value  of  elliptical  circulation  distribution 
= spanwise  coordinate  measured  from  midspan 
3 rotor  radius 


Making  the  transformation  y = R cos  0 , dy  = -R  sin  0 d0,  the  thrust  becomes 


T = PVr^R  sin^  0 d0 


= - PVr  R 
2 o 

Then  the  strength  of  the  trailing  vortices  is  given  by  the  strength  of  the  midspan 
circulation,  i.e.. 


^ PVR 
2 


In  the  for  wake  the  two  vortices  ore  assumed  to  be  convected  downward  at 
one-holf  the  local  downwash  velocity,  due  to  vortex  sheet  rollup  effects,  as  suggested 
by  the  results  on  page  45  of  Reference  34.  The  resultant  flow  velocity  is  then  given  by 


the  vector 


V = ((V  + 2v  cos  i),  (v  sin  i)l 
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where  from  Reference  36,  page  185, 


V 


C^R 


2 l/+X^l 


2,1/2 


and 


U 


V sin  i 

nR 


V cos  i -t-  V 

nR 


(5) 


(6) 

(7) 


c 


T 


T 

pnR2  (nR)2 


(8) 


Consequently,  the  wake  deflection  ongle  x to  the  freestreom  is 


(91 
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INCOMPRESSIBLE  POTENTIAL  AERODYNAMICS  FOR  SEPARATED  FLOWS 

The  general  formulation  for  incompressible  potential  aerodynamics  for 
separated  flows  is  presented  in  this  section.  For  incompressible  fluids,  potential  flows 
are  obtained  if  the  fluid  is  inviscid  and  the  velocity  flow  field  is  irrotational  at  time 
t = 0.  Accordingly,  these  hypotheses  are  invoked  here.  Since  the  analysis  of  sepa- 
rated flows  was  not  included  in  the  preceding  descriptions  of  the  formulation  for 
SOUSSA,  a new  formulation  is  presented  here,  starting  from  basic  principles  and 
avoiding  the  hypothesis  that  the  flow  is  not  separated. 


3.1  INCOMPRESSIBLE  INVISCID  FLOWS 


Assuming  the  fluid  to  be  [ ll*  inviscid,  and  [21  incompressible,  the 
motion  is  governed  by  the  Euler  equations 

[l;  + V - vj  V = - -L  Vp  (10) 

and  the  continuity  equation  (for  incompressible  fluid) 

V • V = 0 (11) 

where  V is  the  velocity  vector  with  respect  to  a prescribed  frame  of  reference,  p is 
the  pressure,  p the  density  (constant  for  incompressible  fluid),  t the  time,  whereas 


V 


“A  o-EA 

' 9x  ' 


(12) 


is  the  del  operator.  Equations  (10)  and  (1 1)  form  a system  of  four  partial  differential 

equations  for  four  unknowns  V , V , V , and  p , 

X y z 


*Numbers  in  brackets  indicate  new  hypotheses  introduced  in  the  formulation. 
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In  order  fo  complete  the  formulation  of  the  problem,  the  boundary  conditions 
at  infinity  on  the  body  and  on  the  wake  must  be  obtained.  Here  it  will  be  assumed  that 
(3l  the  unperturbed  flow  consists  of  uniform  translation,  e.g. , in  the  x-axis  direction. 
Hence,  the  boundary  condition  at  infinity  may  be  written  as 

V = U^T  (for  P atoo)  (13) 

where  is  the  freestream  velocity.  In  particular,  = 0 for  helicopters  in  hover. 

On  the  body  it  is  assumed  that  (4l  the  surface  of  the  body  is  impermeable, 
i.e. , the  normal  components  of  the  velocity  V of  the  fluid,  and  of  the  velocity  Vg  of 
the  body  coincide  at  point  P on  the  surface  ^ g of  the  body: 

(V  - Vg)  >0=0  (for  P on  Tg)  (14) 

where  n is  the  normal  to  Tg  at  P. 

Next,  the  boundary  condition  on  the  wake  for  lifting  bodies  is  considered. 

It  is  assumed  that  [5]  the  body  has  certain  separation  lines  (in  particular  sharp  trailing 
edges)  and  that  the  velocity  flow  field  has  surfaces  of  discontinuity  (wakes)  emanating 
from  these  separation  lines  (or  trailing  edges).  Indicating  by  subscripts  1 and  2 the  two 
sides  of  the  wake,  let  n be  the  outward  normal  on  side  1 and  let 

Af  = 

denote  the  discontinuity  of  any  functions  f acrou  the  wake  surface.  (For  the  classical 
wing  wake,  1 and  2 correspond  to  upper  and  lower  sides,  respectively;  n is  the  upper 
normal  and  Af  = f^  - f^.)  The  boundary  condition  on  the  wake  surface  is 
that  there  Is  no  pressure  discontinuity  at  point  P on  the  woke 

Ap  = 0 (for  P on  E^)  (16) 
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This  condition  con  be  expressed  more  conveniently  in  terms  of  the  velocity  discon- 
tinuity. However,  this  requires  the  introduction  of  the  concept  of  Lamb  surfaces. 
These  surfaces  may  be  avoided  in  the  case  of  potential  flow  and,  therefore,  the 
boundary  condition  is  given  here  for  potential  flows  only. 

Equations  (10),  (11),  (13),  (14),  and  (16)  completely  define  the  problem 
in  terms  of  four  partial  differential  equations  for  four  unknowns  V^,  V^,  and  p, 
with  the  corresponding  boundary  conditions. 

3.2  DIFFERENTIAL  FORMULATION  FOR  INCOMPRESSIBLE 

POTENTIAL  FLOWS 

In  the  preceding  section,  the  flow  was  assumed  to  be  1 1 1 inviscid  and  [ 2] 
incompressible.  Next  assume  ( 6]  the  flow  to  be  initially  irrotational;  then  the  flow 
remains  irrotational  at  all  times  and,  therefore,  there  exists  a potential  function  f 
such  that 

V = V « (17) 

In  this  case.  Equation  (lu/  may  be  integrated  to  yield  Bernoulli's  theorem* 

|i  .i  (18) 


*The  x derivative  af  Equation  (18)  yields 


9x  \ 9t 


L 1 

|9»I^  + 

= - 1 

a*  \ a*  a^* 

2 

at  i 

ax/  ax  dx2 

ijL  + 1 ^ 

az 

axaz  p 9x 

a# 


= (i.  + 9.V  V + 1 iE  = 0 

\at  / p ax 


in  agreement  with  the  first  component  of  Equation  (10). 
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The  constant  on  the  right-hand  side  of  Equation  (18)  is  obtained  from  the  boundary 


conditions  at  infinity. 

Furthermore,  Equation  (11)  may  be  rewritten  as 

= 0 

where  V is  the  Laplacian  operator. 

Similarly,  Equation  (13)  may  be  rewritten  as 

♦ = U X (for  P at  oo) 

00 

and  Equation  (14)  becomes 


(19) 


(20) 


a*  w 

— = Vj  • n 

on 


(21) 


Next  consider  the  boundary  condition  on  the  wake  given  by  Equation  (16). 
Combining  Equations  (16)  and  (18)  yields 


or 


wi 


ith 


and 


9 / . - ' 1 f V« , • 7* , - = 0 


at  ' ' ^’2 


1 ’'1 


Li.  + V*  • =0 

Ut  A 


= » 


1 ’2 


V,  = 1 1 V, 


(22) 


(23) 


(24) 


(25) 


2 \ • ‘/  2 \ ■ -/ 
is  the  average  value  of  the  velocity.  This  average  value  may  be  attributed  to  the  point 
P on  the  wake.  Equation  (23)  may  thus  be  interpreted  os 


w 


Li  (P  ) = constant  in  time  (26) 

NV 


where  is  the  physical  point  of  the  woke  being  convected  with  velocity  V^.  Note 
that  Equation  (23)  implies  that  the  wake  is  tangent  to  V^,  In  particular,  for  the 
steady  state,  the  wake  is  composed  of  the  streamlines  of  points  emanating  from 
trailing  edges  or  other  separation  lines.  Equation  (19)  is  a partial  differential 
equation  for  i with  boundary  conditions  given  by  Equations  (201,  (211,  and  (231. 
Once  $ is  known,  the  velocity  is  given  by  Equation  (17),  while  the  pressure  is  ob- 
tained from  Equation  (18).  It  may  be  worth  noting  that  the  potential  aerodynamics 
formulation  considered  in  this  section  is  considerably  simpler  than  the  formulation 
based  upon  physical  variables  (V  and  p)  since  the  former  gives  one  linear  partial 
differential  equotion  insteud  of  four  nonlinear  ones. 

3.3  PERTURBATION  POTENTIAL 

In  order  to  use  the  Green's  theorem  approach,  it  is  convenient  to  have 
homogeneous  boundary  conditions  at  infinity.  This  can  be  accomplished  by  introduc- 
ing the  perturbation  potential  cp  , related  to  the  potential  f by 

» = U-^  X 4 <p  (271 

Then,  Equation  (19)  yields 

=0  (281 
while  the  condition  at  infinity.  Equation  (201,  becomes  homogeneous 

tp  = 0 (for  P at  CO).  (291 

On  the  other  hand,  combining  Equations  (211  and  (271,  the  boundary 
condition  on  the  body  is  now  given  by 
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(30) 


^ = f-U_  i +Vgj  • n 


dn 


CX) 


while  on  the  woke  (see  Equations  (23)  and  (26)) 


-i-  + V*  -V  6tp  =0 


(31) 


lat 


or 


Acp  (P^)  = constant  in  time 


(32) 


Furthermore,  once  9 is  known,  the  perturbation  velocity 


V = 7«p 


(33) 


may  be  evaluated  and  then  the  velocity  is  given  by 


V = U^i  + v. 


(34) 


Accordingly,  the  Bernoulli  theorem  may  be  rewritten  as 


^ — + - l7«Pl^  + -E  = 


at 


'ao 


dx  2 


Pco 


(35) 


The  above  equations  could  be  obtained  in  a slightly  different  way,  i.e. 
by  assuming  the  frame  of  reference  to  be  connected  with  the  undisturbed  air.  The 
velocity  of  the  body  in  such  a frame  of  reference  would  be  given  by  -U^  ' 
in  agreement  with  Equation  (30). 


3.4  GREEN'S  FORMULA  FOR  LAPLACE  EQUATION 

Consider  two  arbitrary  functions  f and  g and  note  that 
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T1 


and 


V • (fVg)  = Vf  • 7g  + f7^g 

7 • (gV  f)  = 7g  • 7f  + gV^f 

(36) 

7 • (f7g  -g7f)  = fy^g  - gV^f . 

(37) 

Also,  according  fo  the  Gauss  theorem,  for  any  arbitrary  vector  a 


V • a dV 


a • n dr 


E 


(38) 


v/here  E surrounds  V'ond  n is  the  inward  directed  normal . Combining  Equations 
(37)  and  (38)  yields 


III  (f  7^g  - g7^f)  dV  = ^ (f?  g - g 7f)  ' " ^E 


Next,  assume  that  both  f and  g satisfy  the  Laplace  equation,  that  is 
T^f  = T^g  = 0 (in  V). 


(39) 


(40) 


Then  Equation  (39)  yields 


(41) 


which  is  the  desired  Green's  formula  for  the  Laplace  equation. 

Next  choose  the  function  g to  be  the  simplest  function  that  satisfies  the 
Laplace  equation.  If  g is  a function  of  the  distance  r from  a specified  point  P«,  then 


lilMaMMJMiAe 
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the  function  that  satisfies  the  Laplace  equation  is 


where 


r = |P  -P*l 


(42) 


(43) 


Note  that  P*  must  be  outside  the  volume  V;  otherwise  Equation  (40)  is  not  satisfied 
at  P*.  Also  note  that  the  value  of  B does  not  affect  Equation  (41)  under  the  condi- 
tion 


(rt  ^ dr  = 0 

^ an 

which  is  obtained  from  Equation  (41)  with  g = 1.  Therefore,  the  convenient  value 
B = 0 is  used  here.  For  this  case  the  function  g represents  a source.  It  is  convenient 
to  choose  A such  that  g represents  a unit  source,  that  is,  a source  with  flux  equal  to 
1.  This  is  obtained  with  A = -1/4n,  i.e.,* 

9 = --J-  <■"' 

4n  r 

Combining  Equations  (41)  and  (44)  yields  the  desired  Green's  formula  for 
the  Laplace  equation 


(46) 


* For,  in  this  case,  the  flux  through  a spherical  surface  of  radius  R is  given  by 

IT  ii  dr  = • 4ttR^  =1 

dr  p 4ttR2 

r r 
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3.5  GREEN’S  THEOREM  FOR  INCOMPRESSIBLE  POTENTIAL  AERODYNAMICS 

In  order  to  obtain  the  Green's  theorem  for  incompressible  aerodynamics 
(Laplace  equation  for  the  exterior  problem,  i.e.,  for  V outside  the  surface  Eq),  it  is 
convenient  to  consider  the  surface  E composed  of  three  branches  as  indicated  in  Figure 
5.  Note  that  the  inward  normal  for  V is  outwcrd  to  is  a spherical  surface  of 

radius  e and  center  P*,  while  ^2  is  a spherical  surface  of  radius  R and  center  P*. 

Then  Equation  (46),  for  f = cp,  reduces  to 


Eo  + ^l  +^2 


cpl_  i dE  = 0 

3n  \4nr/  dn  \4rtr/ 


Note  that  as  the  radius  e of  goes  to  zero 


n w \4^r/  dn  \4nr  / 

£-•0  E^  L V ' ' 

. ,i„  L (ft  -± 

£-0  L E,  \ ' ' ' 


= k^(ftdr  Jii) 

? ^ \dn/*  4n£  •• 


0 L 4n£^  E 


= lim  9*  + — e = «P* 
« I ^ n , * 


whereas,  as  the  radius  Rr-  of  Eo  goes  to  infinity 

2 
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I 


I 


(49) 


under  the  conditions* 


lim  <p  = 0 
R-^ 

"2 

lim  ^ R = 0 


R-<30 

^2 


dn 


(50) 

(51) 


Therefore,  os  < tends  to  zero  and  tends  to  infinity,  Equation  (37) 


becomes 


(52) 


If  the  point  P*  is  inside  Eq,  then  there  is  no  need  for  the  surface 
In  this  case.  Equation  (47)  yields  simply 


* It  is  verified  o posteriori  (Equations  (71)  and  (72))  that  these  conditions  ore 
satisfied. 

)1 


I 


where 

E*  = E(P*)  = 1 (P*  outside  Eq) 

= 0 (P*  inside  Eq)  (55) 

3.6  ROTOR-FUSELAGE-WAKE  INTERACTION  WITH  FLOW  SEPARATION 

The  model  considered  here  for  the  analysis  of  separated  flows  around  rotor- 
fuselage  configurations  consists  of  the  actual  geometry  of  the  fuselage  and  the  rotor 
blades.  In  addition,  o vortex-layer  wake  emanates  from  the  trailing  edge  of  each 
rotor  blade,  and  a vortex -layer  wake  emanates  from  the  separation  line  os  shown  in 
Figure  6.  Finally,  the  vorticity  in  the  separated  wake  is  modeled  as  one  single 
isolated  vortex  line. 

The  model  described  above  is  relatively  crude.  In  particular,  isolated 
vortices  are  transported  by  the  flow  field  whereas  the  isolated  vortex  considered  here 
is  assumed  to  remain  in  a fixed  position  with  respect  to  the  helicopter.  This  apparent 
contradiction  is  explained  by  the  fact  that  the  distributed  vorticity  in  the  separated 
wake  (converted  downstreom  by  the  flow  field)  is  continuously  replaced  by  the  vorticity 
coming  from  the  boundary  layer.  Therefore,  the  vortex  line  should  be  thought  of  as  an 
analytical  model  to  represent  the  vorticity  field  in  a form  that  is  computationally 
efficient.  As  a consequence,  the  location  of  the  vortex  line  is  somewhat  arbitrary, 
while  the  intensity  can  be  used  to  obtain  the  condition  that  the  velocity  be  equal  to 
zero  at  (Figure  6). 

The  model  described  above  can  be  analyzed  using  the  potential  flow  formulo- 
tion  presented  in  the  preceding  section.  In  order  to  accomplish  this,  note  that  the  flow 
field  is  potential  everywhere  except  on  the  wakes  and  riie  isolated  vortex  line.  Tliere- 
fore,  if  one  considers  the  surface  Eq  depicted  in  Figure  7 (i.e.,  a surface  surrounding 
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the  rotor,  fuselage,  wakes,  and  the  isolated  vortex)  the  flow  is  potential  outside  Zq 
and,  therefore.  Equation  (54)  is  valid. 


Next  consider  an  open  surface  having  as  its  boundary  the  isolated  vortex 
line.  Such  a surface  does  not  intersect  any  of  the  other  surfaces  (i.e. , the  fuselage, 
rotor,  rotor  wake,  and  separation  woke)  and  has  otherwise  arbitrary  shape.  The  func- 
tion of  this  surface  is  similar  to  the  one  of  the  branch  lines  in  the  theory  of  complex 
variables  and,  therefore,  it  will  be  indicated  as  a branch  surface,  or  bronch  wake, 

Eg,  as  shown  in  Figure  8. 

Next  let  Eq  become  infinitesimally  close  to  the  surfaces  of  the  fuselage 
and  wakes  (including  the  branch  wake).  Note  that  in  this  process,  the  closed  surface, 
E^,  surrounding  each  wake  is  replaced  by  the  two  sides  of  an  open  surface  E^  such 
that  E'yy  = E'l^  + E'g.  For  this  case 


(ff  i?.ldE-ff  dE  = 0 


dn  r 


sinceAl— ) =0,  whereas 


Therefore,  in  the  limit.  Equation  (54)  reduces  to* 


4nE*«P*  = -(ft  — 1 -9  — f-i]  dE  + rr  A9  — (i]dE 
^ an  r an  \ r / an  \ r i 


Ef+Er 


E'f+E*b+E-r 


Note  that  the  wake  surfaces  E'p,  E'o,  and  E'|^  are  open  surfaces.  The  source 
integrals  on  the  wake  surfaces  are  identically  equal  to  zero  since  the  values  of 
acp^n  on  the  opposite  sides  of  these  surfaces  are  equal  to  zero. 


i' 

where  Zp  is  fhe  (closed)  surface  of  the  fuselage,  is  the  (clased)  surface  of  the  rotor 
blades,  Z'p  is  the  (open)  surface  of  the  body  wake  (emanating  from  the  separation  line), 

Z'g  is  the  (open)  surface  of  the  branch  wake,  and  Z'j^  represents  the  wakes  of  the  rotor 
blades  (see  Figure  8).  Furthermore  (see  Equation  (15)), 

= 9^  -«P2  (59) 

while  n is  the  normal  on  the  side  1 of  the  wake.  Note  that  A9  is  known  from  Equation 
(32),  on  Z ' and  Z'  , while  A9  is  constant  and  equal  to  the  vortex  intensity  on  Z'„ 

(Appendix  A).  Note  also  that  the  vortex-layer  wake  of  the  fuselage  and  rotor  are 
represented  as  a doublet  layer.  The  proof  of  the  equivalence  of  doublet  loyers  and 
vortex  layers  is  briefly  outlined  in  Reference  40. 

In  addition,  if  P*  is  on  Zp  + Zj^  then,  as  shown  in  Reference  23,  E*  = 1/2. 

Thus,  Equation  (55)  may  be  rewritten  as 

E*  = E(P*)  = 1 (P*  outside  Zp  + Z|^) 

= 1/2  (P^onZp+Z^) 

= 0 (P*  inside  Zp  + Zj^)  (60) 

Note  that  on  Zp  and  Zj^,  ?cp/9n  is  prescribed  by  the  boundary  conditions  of  Equa- 
tion (30).  Thus  for  P*  on  Zp  + Z[^,  Equation  (54)  is  an  integral  equation  relating  cp  to 
the  prescribed  ?9/5n.  That  equation  is  used  to  analyze  the  potential  flow  for  inter- 
action of  fuselage,  rotor,  and  wakes  in  the  presence  of  flow  separation. 

3.7  ALTERNATIVE  TREATMENT  OF  THE  ISOLATED  VORTEX 

An  alternative  treatment  of  the  isolated  vortex  is  presented  here.  This 
treatment  is  similar  ta  the  formulation  of  rotational  aerodynamics  presented  in  Refer- 
ence 30,  and  the  derivatian  of  this  treatment  could  be  obtained  from  that  reference. 
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However,  since  the  vorticity  Is  represented  by  a single  vortex,  the  some  results  may 
be  obtained  directly  from  Equation  (58). 


Note  that  combining  Equations  (33)  and  (58),  one  obtains 


V = 7^9  = 7*(J1)  ^ (^1  - <p  —(—1  dZ 

A -N  j 


dn  \4n  r / dn  \4n  r 


— /I  ] dr 


dn  I r 


ri  j-r'  4-r' 

F ^ B ^ R 


Equation  (61)  may  be  rewritten  as 


V = 7,^  (p  + V. 


where  Vj,  is  the  velocity  induced  by  the  isolated  vortex  (with  = T;  see  Equations 
(A-26)  and  (A-27)  of  Appendix  A), 


V = r7*  IT  A _A]  dr 
‘ JJ  dn  4TTr/ 


while is  the  potential  in  absence  of  the  isolated  vortex,  which  therefore  satisfies 


the  equation 


4,E.i.  is  U 

dn  v4tt  r / dn  v4tt  r 


11  ^ (f) 


r»  +E* 


Nofe  that  the  boundary  condition  on  ^ is  different  from  the  one  on  «P,  for,  according 
to  Equations  (14),  (34)  and  (40),  the  boundary  condition  for  9 is 


(-U^  I - Vj,  + Vg)  . n 


Note  that  in  this  formulation  Vg  appears  only  in  the  boundary  conditions.  Also, 
according  to  Equations  (62),  (A-27),  and  (A-30), 


Equation  (66)  implies  the  evaluation  of  a vortex  line  integral  while  Equation  (63) 
(and  thus  the  preceding  formulation)  requires  the  evaluation  of  a doublet  layer  inte- 
gral. The  last  is  much  more  lengthy  than  the  first  one.  Therefore,  the  formulation 
presented  here  is  preferable  for  the  analysis  of  the  problem  considered  here. 

A further  odvantage  is  that  the  formulation  may  be  easily  extended  to  the 
case  in  which  the  vorticity  is  represented  by  a continuous  distribution.  For  this  case 
this  formulation  is  still  valid  if  Equation  (66)  is  replaced  by  (see  Reference  30) 


Vp  = 7^x 


III  ‘ .ir 


For  the  above  reasons  the  formulation  presented  in  this  subsection  is  used  in 
the  remainder  of  this  report. 

3.8  CONDITIONS  AT  INFINITY 

In  this  subsection  Equations  (50)  and  (51)  are  verified.  Note  that  as  r 
goes  to  infinity.  Equation  (52)  yields  \i 
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SECTION  4 
WAKE  DYNAMICS 


In  this  section,  the  fundamental  aspects  of  wake  dynamics  are  presented. 
These  require  some  basic  concepts  on  the  dynamics  of  vortex  layers.  Also,  it  should 
be  noted  that  vortex  layers  are  not  used  in  Program  SOUSSA,  which  deals  only  with 
source  and  doublet  layers.  However,  it  may  be  shown  that  a vortex  layer  (and  isolated 
vortices)  may  be  represented  by  a doublet  layer.  Therefore,  discussion  on  doublet 
and  vortex  layers  is  presented  in  Appendix  A as  background  material;  source  layers 
are  also  included.  Here  vorticity  dynamics  are  presented  and  finally  the  results  are 
applied  to  the  wake  dynamics. 


4.1  VORTICITY  DYNAMICS 

In  order  to  formulate  the  wake  dynamics  problem,  it  is  necessary  to 
reformulate  some  classical  results  on  vorticity  dynamics  in  a novel  fashion  that  is 
particularly  suitable  for  the  analysis  of  wake  dynamics. 

Consider  the  equilibrium  equation  for  an  inviscid  fluid  (Euler's  equation) 

DV  _ _ 1 VP  r 


Note  that 


DV  ^ 

Di  at 


— + 1 r (V-V)  + C X V 
St  2 
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where 


C = ^ X V 


(75) 


Taking  the  curl  of  Equation  (73)  yields,  for  incompressible  flow 


f X 2V  = 0 

Dt 


or,  according  to  Equations  (74)  and  (75) 


— + Vx(CxV)  = 0 

at 


(76) 


Equation  (76)  may  be  rewritten  in  a more  useful  form  by  noting  that 
(Reference  41 , page  43) 

? X (axb^)  = (b*V)d'  - (a-v)b  - b(r-a)  + a(r-b) 


(77) 


and,  therefore,  combining  Equations  (76)  and  (77)  and  noting  that  V • C 0 ^ 

= 0 (incompressible  flow),  yields 


— + (V-V)C  - (C-V)V  ^ 0 

at 


(78) 


or 


— = (C-r)V 
Dt 


(79) 


Equation  (79)  is  a classical  result  known  as  the  third  vortex  theorem. 
Equation  (79)  may  be  integrated  in  closed  form  as  shown  in  Reference  41 , page  80. 
However,  a novel  method  of  solution  was  obtained  in  Reference  30  and  is  briefly 
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presented  here  since  this  new  result  facilitates  the  discussion  on  wake  dynamics. 
Note  first  that  by  introducing  a system  of  curvilinear  coordinates,  , the  vector  C 
may  be  written  as  (Reference  30) 


- Of  — 

C ~ C 9a 


(80) 


where 


ax 


a = 1,  2,  3 


(81) 


are  the  base  vectors,  while  are  the  contravariant  components  of  the  vector  Q.  Note 


that 


= C-g" 


(82) 


where 


-a  _ - 

g - 9 gp 


(83) 


Wl 


ith 


orp  _ o 

g gpY  °Y 


(84) 


Of 


where  6 ^ is  the  Kronecker  delta  and  is  the  metric  tensor  given  by 


gpY 


(85) 


Note  that 


Dt 


^ 9.  - 


Dt 


Dt 


(86) 
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Nexf  assume  that  the  coordinate  system,  moves  with  the  fluid.  Such  a system 
describes  the  motion  from  a Lograngian  (rather  than  Eulerian)  point  of  view;  therefore, 
it  may  be  called  a Lograngian  coordinate  system.  Such  a system  is  also  called  a materi- 
al  coordinate  system  (since  at  any  time  the  point,  § = constant,  represents  a material 
point)  or  convected  coordinate  system  (since  the  coordinate  points  are  convected  by  the 
system).  Note  that,  in  this  coordinate  system,  the  total  time  derivative  coincides  with 
the  partial  time  derivative,  since  the  total  time  derivative  is,  by  definition,  the  time 
derivative  in  a system  moving  with  the  fluid.  Tnerefore, 


(87) 


(88) 


(89) 


(90) 


I .e. , 


c“  (5’,  t)  = c“(?^  0)  (92) 

O' 

Equation  (92)  indicates  that  in  a Lagrangian  coordinate  system,  Q is 
ot  a 

only  a function  of  § ; i.e.,  C 's  a function  of  the  material  point  and  not  of  time. 

a 

In  particular,  in  steady  state,  C is  constant  along  a streamline. 

The  result  obtained  above  yields  two  important  properties  of  flow  fields: 

(a)  A particle  or  a portion  of  fluid  initially  irrotational  remains 
irrotational  at  all  times  (Lagrange-Cauchy  theorem). 

(b)  Vortex  lines  are  material  lines^  i.e.,  the  set  of  particles  that 
composes  a vortex  line  at  one  instant  will  continue  to  form  a vortex 

line  at  later  instants. 

The  first  property  is  an  immediate  consequence  of  Equation  (91),  In 
order  to  prove  the  second  one,  consider  a material  line,  L^,  described  by  a parameter 

"X  = X (?"(H),  tl  = X(>v  t)  (93) 

which  at  time  t = 0 coincides  with  a vortex  line,  i.e.,  C 's  tangent  to  L^,  or 

C(H,  0)  = 01  = \ {»<)^|?“(H),  0)  (94) 

Then  the  proof  of  property  (b)  consists  of  showing  that  C remoins  porallel  to  at 
all  times.  Note  first  that  Eauation  (94)  is  equivalent  to 

C“l  !"■(«),  0|g;,|5“(«),  0|=X  ^ if. 

C?  OH 

= \q(m)^  (M)  ^|P/(H),  0)  (95) 

3 H 


c“[5“(h),  0]  = \q(>^  (H) 

dn 


At  time  t,  using  Equations  (92)  and  (96)  results  in 


= C“(?“(H),  0]  i;,[l“(H),t] 


an 

= ^ [?"(>^)/  n 


Equation  (97)  indicates  that  J remains  parallel  to  the  material  line,  L^, 
and,  therefore,  vortex  lines  are  material  lines. 

An  additional  result  (on  vortex  stretching)  may  be  obtained  from  Equation 
(97).  Consider  an  element,  d^',  of  the  vortex  line,  given  by 

dl  = ^ dH  (9£ 


Then  Equation  (97)  indicates  that 

I C I = I C lo 

IdAl  |d£|^ 

which  is  the  well-known  vortex-stretching  law  for  incompressible  fluid;  the  ratio 


I 
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between  the  Intensity,  I ^ | , of  the  vorticity  and  the  length,  jdj^l , of  the  vortex  element 


remains  constant  in  time/ 


WAKE  DYNAMICS 


In  the  formulation  presented  in  Section  3,  the  wakes  are  represented  as 
doublet  layers.  On  the  other  hand,  it  is  shown  in  Appendix  A that  doublet  layers  are 
equivalent  to  vortex  sheets.  Also,  a vortex  sheet  may  be  thought  of  as  the  limiting 
use  of  a distribution  of  vorticity  in  a thin  volume  surrounding  the  vortex  sheet.  This 
is  very  close  to  the  physical  reality;  the  wake  is  actually  the  volume  of  the  flow 
field  where  the  vorticity  generated  by  the  presence  of  the  aircraft  is  not  negligible. 

The  above  considerations  indicate  that  the  results  obtained  in  Subsection 
4.1  are  applicable  to  the  analysis  of  the  dynamics  of  the  doublet  layer  wake.  The 
first  consequence  is  that  the  points  of  the  wake  may  be  considered  as  material  points. 
However,  the  velocity  flow  field  is  discontinuous  through  a doublet  layer  and,  there- 
fore, the  velocity  of  a point  of  the  wake  is  not  clearly  defined.  It  was  Helmholtz  who 
first  suggested  that  the  velocity  of  a point  on  the  wake  is  the  average  of  the  values  on 
the  two  sides  of  the  wake 


A 2 


This  result  may  be  obtained  by  evaluating  the  velocity  of  a point  on  the  wake  in  the 
limiting  use  as  the  thin  layer  of  vorticity  reduces  to  a vortex  sheet.  Also,  this  is  in 
agreement  with  the  result  obtained  in  Section  3 where  it  was  shown  that  (see  Equation 

(31)) 

*Note  that  the  same  result  is  obtained  from  Kelvin's  theorem,  which  yields  in  the  limit, 
r = ^ 7-.d5  = QdV/dJi  ~ constant 
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D^(&cp ) 


0 


(101) 


where 


may  now  be  interpreted  as  the  total  derivative  for  a material  point  on  the  wake. 

Based  on  the  above  discussion,  the  following  procedure  may  be  used  to 
evaluate  the  geometry  of  the  wake.  First,  a point  on  the  wake  is  a material  point 
having  velocity,  V^.  Therefore,  given  its  location  at  this  time,  t^,  the  new 
location  at  time  t + dt  is  given  by 


This  yields  the  geometry  of  the  wake  at  any  time  t.  Furthermore,  Equation  (101) 
indicates  that  A'f  is  a function  of  but  not  of  time.  Therefore,  the  distribution  of 
iP  on  the  wake  is  immediately  available  once  the  location  of  each  material  point  is 
obtained. 

Summarizing,  the  dynamics  of  the  wake  is  solved  by  following  the  loca- 
tions of  the  material  points  in  time  and  by  the  fact  that  using  the  value  of  A'P 
corresponding  to  each  material  point  does  not  change  in  time. 


SECTION  5 


NUMERICAL  FORMULATION 


The  numerical  formulation  used  to  approximate  Equation  (58)  with  a system 
of  algebraic  equations  is  presented  in  this  section.  In  summary,  the  procedure  is  as 
follows.  The  surfaces  of  the  rotor,  fuselage  and  their  respective  wakes  are  divided  into 
small  surface  elements,  E..  The  potential  V and  the  normal  wash  V are  assumed  to  be 
constant  within  each  element.  This  yields  a set  of  linear  delay  equations  relating  the 
values  of  the  potential,  at  the  centroids  of  the  elements,  E^,  to  the  values  of  the 
normal  wash,  at  the  centroids  of  the  elements,  E^.  The  coefficients  are  evaluated 
analytically  for  quadrilateral  hyperboloidal  elements. 

5.1  APPROXIMATION  OF  INTEGRAL  EQUATION 


For  the  geometry  used  in  deriving  Equation  (58)  (see  Figure  8),  divide 

the  surfaces  of  fuselage  and  rotar  Ej.  and  E„,  into  N surface  elements  E . Also,  divide 

I K n 

the  surfaces  of  their  wakes,  E,.'  and  En'  into  M surface  elements E '.  Then  Equation 

r K m 

(58)  may  be  rewritten  as: 


2E(P*,t)cp(P*,t) 


dE 


n 


(104) 
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2E(P*/t)'P  (P*»0  = 

n = 1 


+ 53  Cn  (P*  ,t)  <Pn  (0 

n = 1 


+ 53  F_  (P*,t)^Vn,  (t) 


m = 1 


? CP 

* = — ^ 

" 9 0 


P = Pn 


8„(P.,0  =0 


En(t)  2TTr 


By  imposing  the  condition  thot  Equotion  (105)  be  satisfied  at  the  centroids 
P,  of  the  elements one  obtains  (note  that  according  to  Equation  (60)  E{Pj,)  - 1/2, 


since  P^  is  on  Zp  + Zj^), 


N N M 

*1,"'  ' Z “hnC' * Z W-'-PnC'^Z 

n=l  n=l  m=1 


where 


Cu„(0  = C„  (P.  ,t) 


It  should  be  emphasized  that  (unlike  in  the  formulation  for  aircraft.  Refer- 
ence 42)  the  coefficients  / ‘^nd  are  time  dependent.  This  is  due  to  the 

fact  that  the  surfaces  andS^  vary  with  time  (except  for  the  ones  of  the  fuselage) 
and  also  to  the  fact  that  the  points  move  in  time  (except  for  those  on  the  fuselage). 

Next  the  relationship  between  Acp  ond  «P  is  obtained.  Consider  the 

m n 

trajectory  of  the  centroid,  P'  , of  the  element z'  of  the  wake.  Let  P^^  indicate  the 

mm  m 

point  of  the  trailing  edge  from  which  the  above  trajectory  starts  (see  Figure  9). 

Then,  according  to  Equation  (32), 


LV  (t)  = ACP  (t  - T ) 
m'  ' m ' m' 


Figure  9.  Relation  of  Wake  Trajectory  to  Trailing  Edge. 
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where  is  the  value  of  Atp  at  whereas  t is  the  time  necessary  for  the 

m mm 

material  point  of  the  wake  to  move  from  to  P^  with  velocity  given  by 
Equation  (25). 

In  addition,  it  is  possible  to  express  the  value  of  'ri  terms  of  the 

values,  9 , of  the  potentiol  at  the  centroids  of  the  elements,  A simple  expression 

. TE 

is  obtained  by  assuming  that  (in  view  of  the  Kutta  condition)  the  volue  of  may 

be  approximated  with  the  value  of  AV  at  the  centroids  of  the  elements  in  contact  with 
the  point  pj^^.  This  yields 

where  S is  equal  to  1 (-1)  for  the  point  of  the  upper-surface  (lower -surface)  element 
in  contact  with  the  point  P^^  of  the  trailing  edge  and  equal  to  zero  otherwise.  In 


other  words, 


S = 1 if  pj^^  (on  upper  surface) 


= -1  if  P €E  (on  lower  surface) 
m n 


= 0 otherwise 


Combining  Equations  (108),  (110),  and  (111)  yields 


N N 

'Pn(0  =E  + E C hn  (^)  'Pn 

n = 1 n = 1 


M N 


L.  E (t-j 


m = 1 n = 1 
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This  is  a system  of  linear  delay  equations  with  time-dependent  coefficients 
that  may  be  solved  with  a step-by-step  procedure,  once  the  coefficients 
C^^(t),  and  have  been  evaluated.  An  explicit  expression  for  these  coefficients 

is  given  in  Subsection  5.2.  An  alternative  iterative  procedure  for  the  solution  of 
Equation  (113)  is  presented  in  Subsection  5.3, 


i 


5.2  ANALYTICAL  EXPRESSION  FOR  COEFFICIENTS  B,  , C,  , F. 

hn  nn  hm 

Consider  the  coefficients  ^hn^^^  that,  according  to  Equation 

(107)  and  (109),  are  given  by 


(114) 


r • n dE 

n 


(1151 


where  r = | r | with 

r - P(t)  -P^(t).  (116) 

P^(t)  is  the  centroid  of  the  element  E^(t),  whereas  P(t)  is  the  dummy  point  of  integra- 
tion on  the  element  E^(t).  The  integrals  in  Equations  (1 14)  and  (1 151  may  be 
evaluated  analytically  if  the  surface  element  T;^(t)  is  quadrilateral  and  is  approximated 
by  the  hyperboloidal  element  of  Reference  43,  i.e.,  the  portion  of  a hyperboloidal 
paraboloid  of  the  type  (see  Figure  10), 

RO  = Pq(0  + ?p^(t)  + T|p^(t)  + ?:  T|p2(f)  (117) 
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where  the  surface  coordinates  ^ and  T]  are  bounded  by 


■1  < < 1 


■1  - T1  - \ 


The  vectors  Pq(0,  P|(0,  92^^)  > P3(^^  °''®  '"elated  to  the  corner  points 

Ppp^^^'  Pmp^'^'  element  E^(t)  by  the  relationship 


1 1 


1 -1 


-1  1 


-1  -1 


Pr,_(0 

PP 

P (0 

^pm' 

P (0 

^mp' 

P (0 


Note  that 


d = a^  X 02  I d ? d T] 


where  a.  (i  = 1 , 2)  are  the  surface  base  vectors 


a,  - :^P/d? 


02  = .>P/3T1 


Note  also  that  the  surface  unit  normal  n may  be  expressed  as 


0 = (a,  K aj)y  |o,  X Oj 


Combining  Equations  (1 14),  (115),  (118),  (120),  and  (122)  yields 


>'->  = 1 1 


— I a^  X 02  1 d ?;  d T] 
2n  r 
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Note  that  both  integrals  are  of  the  type 


r • Oj  X 02  d T| 


.1  ] 

I I I M?’,  ?:^)dF  dTl  (12 

-1  -1 

which  may  be  evaluated  as 

I F(l,  1)  - F(l,  -1)  - F(-l,  1)  + F(-l,  -1)  (12 

where  F(F , T))  is  such  that 
^2 

(12 

Therefore,  Bj  (t)  and  Ci  (t)  are  easily  evaluated  if  one  notes  that  (Reference  43) 


— r • a,  X a,,  = 

r3  1 2 


^ r X a^ 

• r X 

1 c 

\\r  • 

Oj  X 

1 c 

whereas,  assuming  n to  be  approximately  constant  within 


1 1 — — — — — 1 -1  / ''  ' ^1 
- 1 a.  X a,j  I = -r  X a,  • n — sinh 

^ I ' |oJ  \17xo,| 


- — - 1 -1  I ' ^2 

t r X 02  • n — — sinh  - — _ 

102!  \l~xa2 


+ r . n tan 


^ r X • r X 02 
\1  r |7  • Oj  X 02 


‘This  is  consistent  with  the  hypothesis  that  cp  and  )|i  are  constant  within  T 


r I 


Next  note  that  the  definition  of  the  coefficients  F,  (t)  is  identical  to  the 

hm 

one  for  the  coefficients  Ci  (t)  with  the  exception  tliot  the  snrfoce  elements  v'  ore  on 

nn 

the  wake.  Therefore,  if  the  wake  elements  are  approximated  hy  quadrilateral  hyper- 
Ix^loidal  elements,  the  coefficients  ore  evaluoted  hy  usinq  the  some  pioceduie 
outlined  for  the  coefficients  C, 
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ITERATIVE  PROCEDURE  FOR  SOLUTION  OF  DELAY  EQUATION 


In  order  to  solve  Equation  (1  13\  a step-Liy-step  procedure  may  lie  used; 
however,  this  is  expected  to  be  time  consuming  (because  of  limitations  on  At  due  to 
stobillty  constraints  usually  connected  with  step-by-step  proceduiesV  Therefore,  for 
the  steady  state  (periodic)  solution,  an  alternote  scheme  may  be  usefully  employed. 
This  scheme  is  based  upon  the  following  iterotive  procedure:  for  the  zeio*^'  iteration, 
assume  that  the  variation  of  Acp  in  time  is  small  enough  so  that 


cp(0)  (t  - T ) ~ 9(0^  (t) 


Introducing  this  approximation.  Equation  (113)  reduces  to 


(1301 


[ \n  - S.n("  - f|™(" 

(1311 

This  is  a system  of  algebraic  equations  (with  time  appearing  os  a pnromettMl 
that  may  be  solved  for  ot  discrete  angular  positions,  \ .,  of  the  rotor,  corres- 

ponding to  time,  t.  Xj/n  (where  Q is  the  angular  velocity  of  the  rotorl.  For  the 
successive  iterations.  Equation  (113)  may  be  rewritten  as 


{.(^>(t)} 


)*n<’'t 


( 


S i (t  - • l| 

mn  [In  * 


:\3?) 
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If  the  iteration  scheme  converges,  the  converged  values  of  tp^^^(t)  ore  the  solutions 


for  Equation  (113). 
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EVALUATION  OF  PERTURBATION  VELOCITY  AND  PRESSURE 


Equation  (1 13)  is  o system  of  linear  delay  equations  with  time -dependent 
coefficients.  After  the  potential  at  the  centroid  of  each  element  is  obtained,  the 
perturbatkon  velocities  and  pressure  can  be  evaluated.  As  shown  in  Subsection  4.1, 
the  velocity  can  be  decomposed  as  follows: 


— —a 

V = V a 

ly 


= 1,2,3 


(133) 


where  a denotes  the  adjoint  base  vector  and  v^  is  the  covariant  component  of 


the  vector  v.  Note  that 


V = 
O' 


(134) 


with  representing  the  curvilinear  coordinates.  If  ? and  T]  are  the  surface  curvilinear 
coordinates  and  a^  and  O2  are  the  base  vectors.  Equation  (133)  becomes 


-1^-2^  -3 

V = v^  a + V2  o + a 


(135) 


where 


-1  °2  ^ °3 

a = 


(136) 


°1  ' °3 


-2  °3  °1 

a = — 


and 


a,  . a2  X 03 


(137) 


-3  °1  °2 


°]  °3 


(138) 


S8 


i, 

fa 


where  i,  j,  and  k are  the  unit  base  vectors. 

Finally,  Bernoulli's  equation  for  incompressible  potential  flow,  given  by 
Equation  (18),  is  used  to  obtain  the  pressure,  that  is. 


5.5  EVALUATION  OF  THE  AERODYNAMIC  COEFFICIENTS 

After  the  pressures  have  been  obtained  at  the  centroids  of  the  aerodynamic 
panels,  the  aerodynamic  coefficients  can  be  obtained  by  summing  the  contribution  of 
each  surface  panel.  Hence,  the  lift  on  the  helicopter  is  given  by  the  expression 


= - JJ 

s 

= - I ^ I ^ Cp  a^  xa2  • k d?dy 


(147) 


since 


d?  d 1) 


4 


(148) 


Similarly,  the  induced  drag  is  given  by 


C 


D 


C n • i dE 
P 


= -4 


E 


C (a,  X a, 

Ph  ’ ^ 


(149) 


Note  that  C denotes  the  pressure  coefficient  at  the  centroid  of  the  h element 

Ph 
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RESULTS  AND  CONCLUSIONS 

The  formulation  presented  in  Sections  3 and  4 for  incompressible  potential 
aerodynamics  with  separated  flow  was  incorporated  into  a computer  program  SHAPES  to 
study  the  rotor  wake  effects  on  hub  pylon/flow  separation.  The  program  SHAPES  is  an 
acronym  for  Subsonic  Helicopter  Aerodynamic  Pogrom  with  Effects  of  Separation. 

The  effects  of  the  presence  of  the  rotor  in  its  various  operating  conditions  os  well  as 
the  presence  of  separation  on  the  pressure  of  the  helicopter  components  bereath  the 
rotor  were  investigated.  Extensive  numerical  results  showing  the  flexibility  and 
accuracy  of  the  method  were  obtained  from  the  program  SHAPES.  Also,  comparison 
with  several  existing  results  was  made.  Concluding  remarks  are  given  in  Subsection 

6.3. 

1 

[ 6.1  NUMERICAL  RESULTS 

Convergence  analysis  was  performed  on  an  ellipsoidal  helicopter  configura- 
tion. The  case  of  a modified  BO-105  helicopter  fuselage  configuration  without  rotor 
was  analyzed  to  study  the  effects  of  separation.  Results  for  a single-bladed  rotor 
and  four-bloded  P<H-51A)  rotor  were  obtained  and  compared  with  available  theoretical/ 
experimental  results.  A fuselage  configuration  (Model  1)  with  a prescribed  separation 
wake  and  a single  two-bladed  rotor  with  classical  helicoidal  wake  geometry  were 
treated  separately  and  combined  for  analysis  at  forward  flight,  and  compared  with 
existing  experimental  results. 

6.1.1  CONVERGENCE 

Convergence  studies  were  conducted  to  identify  the  approximate  number  of 
aerodynamics  elements  required  for  the  analysis. 


Convergence  analyses  were  performed  on  on  ellipsoidal  fuselage  configura- 
tion with  separation  effects  at  an  airspeed  of  150  mph,  zero  angle  of  attack,  and  zero 
sideslip  angle.  The  overall  length,  width,  and  height  of  the  fuselage  were  41,  14.5, 
and  16.5  inches,  respectively.  Results  are  given  in  Figure  11.  Two  different  aero- 
dynamic paneling  schemes  were  used:  nonuniform  and  unifarm  elements. 

Results  for  nonuniform  elements  are  given  in  Figures  11(a)  through  11(d); 
Figures  11(e)  through  11(h)  present  corresponding  results  using  uniform  elements.  Both 
the  velocity  potential  and  pressure  distribution  al^ng  the  top  fuselage  centerline  and 
the  side  of  the  fuselage  were  obtained.  The  four  cas5*s  used  for  each  aerodynamic 
panel  breakup  scheme  were  as  follows: 


Run 



Symbol 

Number  of  Elements 

1 

A 

64 

2 

O 

120 

3 

□ 

168 

4 

o 

200 

It  is  noted  that  the  geometric  symmetry  of  the  fuselage,  i.e.,  left-  and  right-hand 
side,  was  not  used,  and  the  number  of  aerodynamic  panels  represents  the  total  number 
of  elements  used.  On  the  CDC  6600  computer  systems,  250  seconds  of  computer  time 
were  required  to  run  four  cases. 

It  is  seen  that  the  use  of  the  uniform  elements  yields  slower  convergence  of 
the  perturbation  potential  near  the  nose  and  tail  sections  of  the  fuselage.  A slower 
rate  af  convergence  is  seen  when  the  pressure  distribution  along  the  top  centerline  and 
the  side  of  the  fuselage  is  computed.  Hence,  the  use  of  nonuniform  elements  is  desir- 
able, i.e.,  smaller  panels  near  the  nose  and  tail  sections  of  the  fuseloge.  It  is 
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Figure  lUc).  Pressure  Coefficient,  Top  Centerline 
(Nonuniform  Elements) 


AXIAL  distance,  INCHES 


PROGRAM  SHAPES 


estimated  that  for  this  fuselage  approximately  168  nonuniform  elements  are  sufficient 
for  convergence,  whereas  the  use  of  the  uniform  element  scheme  requires  more  than 
200  elements  for  analysis. 

6.1.2  SEPARATION  EFFECTS 

The  study  of  the  effects  of  the  separation  wake  on  an  ellipsoidal  helicopter 
fuselage  configuration  is  given  in  Figure  12.  The  location  of  the  separation  line  is 
prescribed  and  can  be  determined  from  existing  flow  visualization  data  or  estimated 
by  the  use  of  analytical  criteria  for  separation. 

The  strength  of  the  isolated  vortex  element  is  voried  such  that  the  velocity 
is  zero  at  the  prescribed  separation  location.  Figure  12  shows  the  effect  of  the 
strength  of  an  Isolated  vortex  element  on  the  potential  difference  distribution  ond 
pressure  on  the  top  centerline  of  a simple  ellipsoidal  fuselage  at  150  mph,  a 0°, 
and  P = 0°.  Again  the  overall  dimensions  of  the  fuseloge  are  41 , 14.5,  and  16.5 
inches,  respectively.  The  four  cases  run  ore  shown  below. 


Run 

Symbol 

Vortex  Strength 

1 

A 

20 

2 

O 

40 

3 

n 

60 

4 

O 

80 

The  number  of  aerodynamic  elements  used  in  each  case  is  186.  On  the  CDC  6600 
computer  system,  the  required  computer  time  was  approximately  80  seconds  per  case, 
i.e.,  for  each  value  af  vortex  strength.  Again,  the  symmetry  of  the  helicopter 
fuselage  was  not  employed. 
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6.1.3  ROTOR  EFFECTS 

Results  for  a single -b laded  rotor  ore  presented.  The  geometry  and  flight 
conditions  of  the  rotor  were  R = 17.5  ft,  C = 1 .083  ft,  Y =2.33  ft,  = 10.61 
with  C-i"  = .00186  and  Q = 355  rpm.  Further  details  are  given  in  Reference  44.  The 
classical  helicoidal  wake  model  is  used  with  the  induced  rotor  velocity  (associated 
with  the  convection  of  the  rotor  wake)  given  by 

The  spanwise  potential  difference  and  the  sectional  lift  are  shown  in  Figure  13.  It  is 
seen  that  90  elements  for  this  configuration  ore  not  sufficient  for  complete  conver- 
gence, although  convergence  is  close.  Reference  45  indicates  that  about  100  panels 
should  be  required  for  analyzing  a single  rotor  blade.  The  cases  considered  are  sum- 
marized as  follows; 


Run 

Symbol 

Number  of  Elements 

1 

o 

42 

2 

o 

64 

3 

□ 

90 

For  each  of  the  above  runs,  the  number  of  wake  spirals,  N^,  used  is  3 and 
the  number  of  surface  elements,  N^,  along  the  flow  direction  is  18.  Further  details 
on  the  convergence  analysis  of  a rotor  configuration  are  given  in  Reference  45.  The 
computer  time  required  is  approximately  400  seconds  for  the  three  cases.  The  results 
were  compared  with  other  theoretical  results,  i.e.,  Landgrebe's  wake  model  and 
classical  wake  model.  It  is  seen  that  the  present  results  are  in  very  good  agreement 
since  the  present  wake  model  and  the  classical  wake  model  do  not  take  into  account 
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spanwise  sectional  lift  Distribution,  lb/ft 
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the  concentration  of  the  vortex  sheet  into  a tip  vortex  (not  included  in  the  present 
analysis).  It  is  noted  that  the  Landgrebe  wake  model  (Reference  46)  is  obtained  from 
experiments  and  includes  the  effect  of  the  tip  vortex.  Details  of  these  results  ore 
presented  in  Reference  44, 

Results  were  obtained  for  the  XH-51A  Helicopter  rotor  at  hand  with  u=  355 
rpm,  <^0=  10.61°,  and  thrust  coefficient  C-j.  = .00574.  The  rotor  radius  is  17.5  ft  with 
root  cutout  at  2.33  ft  and  constant  blade  chord  of  1 .083  ft.  Using  a classical  heli- 
coidol  wake  model,  the  induced  rotor  velocity  is  estimated  to  be  35  ft/sec.  The  four 
coses  considered  are ; 


Run 

Symbol 

Number  of  Elements 

1 

A 

32 

2 

O 

168 

3 

o 

192 

4 

□ 

200 

In  each  cose,  the  values  = 3 and  = 18  were  prescribed.  Note  that  the  number 
of  elements  represents  the  total  number  of  panels  on  the  four  rotor  blades. 

Convergence  analysis  for  this  rotor  configuration  is  shown  in  Figures  14(a) 
through  14(c),  It  is  seen  from  these  studies  that  the  use  of  200  elements  is  very  nearly 
sufficient  for  convergence.  Both  the  spanwise  potential  difference  and  sectional  lift 
distribution  are  obtained.  These  results  are  compared  with  other  existing  theoretical 
and  experimental  results  and  are  shown  in  Figure  14(c).  Results  are  in  excellent 
agreement  with  experimental  data  for  the  XH-51A  obtained  in  Reference  47.  Also, 
present  results  are  in  good  agreement  with  results  obtained  using  both  the  Landgrebe 
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wake  model  and  classical  wake  model.  Overall, the  results  ore  in  exceptionally  good 
agreement  If  one  considers  that  the  effect  of  the  concentration  of  the  vorticity  into  a tip 
vortex  is  not  taken  into  account.  This  effect  is  clearly  shown  by  Landgrebe  wake 
results  since  the  Landgrebe  wake  model  is  obtained  from  experiments  and  includes  the 
effect  of  the  tip  vortex.  About  1000  seconds  were  required  to  run  o single  cose  (200 
elements)  on  a CDC  6600  computer  system. 

In  both  cases,  the  rotor  wake  geometry  is  helicoidal  and  automatically 
adjusts  itself  for  various  values  of  and  where  is  the  induced  rotor  velocity 
involved  with  the  transport  of  the  rotor  wake  and  S?is  the  rotor  angular  velocity. 

6.2  COMPARISON  WITH  EXISTING  RESULTS 

Results  for  the  BO -105  configuration  without  pylon  and  tail  boom  section  were 
obtained  and  compared  with  available  experimental  results.  Details  on  the  experi- 
mental data  are  given  in  Reference  5.  A sketch  of  the  BO- 105  fuselage  section  is 
given  in  Figure  15.  Oil  flow  visualization  on  the  BO-105  configuration  (Figure  16) 
at  150  MPH  and  rv  = 0 was  used  to  estimate  where  separation  occurs.  The  strength  of 
the  isolated  vortex  element  was  varied  until  close  agreement  with  experimental  results 

was  obtained.  The  results  of  this  study  are  given  in  Figure  17.  In  obtaining  most 
of  these  results,  no  attempt  was  mode  to  verify  the  location  of  the  separation  line. 

When  separation  occurs,  a step  pressure  gradient  is  observed  and  the  velocity  changes 
sign.  It  should  be  noted  that  the  location  of  maximum  total  potential  corresponds  to  the 
location  at  which  the  velocity  changes  sign.  This  is  the  criterion  used  to  determine  the 
correct  vortex  strength  for  the  flow  and  the  result  is  presented  in  Figure  17.  In 
general , the  results  are  good  although  some  improvement  is  desirable.  Improvement  con 
he  accomplished  by  moving  the  separation  line,  varying  the  vortex  strength  at  smaller 
• ■•ments,  and  including  a more  complete  vortex  dynamic  model. 
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Next,  the  results  for  Model  1 with  separation  obtained  from  Program  SHAPES 


ore  compared  with  available  experimental  results.  The  sketch  of  Model  1 and  the  rotor 
blade  plonform  ore  shown  in  Figure  18.  Further  details  of  this  model  ore  given  in  Ref- 
erence 49.  The  cose  of  a fuselage  configuration  without  rotor  geometry  was  investi- 
gated. A total  of  64  elements  was  used  with  o required  computer  time  of  31  seconds  on 
CDC  6600  computer  systems.  It  should  be  noted  that  experimental  data  obtained  by 
Wilson  and  Mineck  (Reference  49)  include  the  presence  of  the  pylon,  whereas  results 
in  Figures  19(a)  and  19(b)  are  presented  without  the  pylon.  Pressure  distribution  on  the 
top  centerline  and  bottom  centerline  of  the  fuselage  was  obtained.  Overall,  the  re- 
sults are  in  good  agreement  since  the  present  result  does  not  include  the  presence  of  the 
pylon  section. 

Results  for  Model  1 with  the  pylon  and  rotor  are  given  in  Figure  19(c). 

The  rotor  radius  is  20  inches  with  a root  cutout  of  8 inches.  The  blade  chord  is  3.5 
inches  with  a thickness  ratio  of  0.12.  In  this  analysis  the  effect  of  separation  was 
included.  The  exact  geometry  of  Model  1 designated  in  Reference  49  is  not  used. 

With  a geometric  preprocessor,  the  geometry  is  approximated  by  simple  overall 
geometric  configuration,  i.e.,  ellipsoidal  and  conical  sections.  Pressure  distribution 
along  the  fuselage  top  centerline  and  bottom  centerline  is  presented.  The  operating 
condition  for  this  run  is  30  knots,  o = 0,  P = 0,  and  rotor  angular  speed  of  3600  rpm. 


The  cases  considered  are  summarized  as  follows. 


Run 

Symbol 

Remarks 

1 

□ 

Fuselage  with  rotor  blade;  no  pylon. 

2 

o 

Fuselage  with  pylon  and  rotor  blade. 

3 

A 

Fuselage  with  pylon,  rotor  shaft,  hub,  shank,  and  blade. 

4 

Fuselage  with  pylon,  rotor  shaft,  hub,  and  shank;  no 
rotor  blade. 

Ll  ~ 
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O EXPERIMENTAL  (WILSON  a MINECK) 
+ PROGRAM  SHAPES  (NO  PYLON) 
NUMBER  OF  ELEMENTS  = 64 


Figure  19(b).  Pressure  Coefficient,  Bottom  Centerline 


Figure  19(c).  Pressure  Coefficient,  Top  Centerl 


A typical  computer  time  required  for  o cose  of  Model  1 with  rotor  blade  is  approximately 
600  seconds  on  the  CDC  6600  computer  system  with  184  elements  used.  Again,  the  re- 
sults obtained  here  compare  very  well  with  experimental  data  presented  in  Reference  49. 
It  is  seen  that  the  effect  of  the  pylon  is  small  because  the  actual  pylon  geometry  Is 


replaced  with  a streamlined  ellipsoidal  pylon  section.  Hence,  the  pressure  change  on 

this  fuselage/pylon  section  is  smooth  and  not  abrupt.  , 

6.3  CONCLUDING  REMARKS  ■ 

The  theoretical  formulation  and  corresponding  numerical  procedure  for  the 
study  of  the  total  effect  of  the  presence  of  the  rotor  in  its  various  operating  conditions 

on  the  pressure  distribution  and  drag  (induced)  of  the  helicopter  components  beneath  the  | 

rotor  were  presented.  A potential  flow  aerodynamic  program,  SHAPES,  with  suitable  rotor  j 

wake  representation  was  developed  to  predict  the  separation  characteristics  of  arbitrary  | 

three-dimensional  helicopter  configurations,  with  lifting  surfaces  in  yawed  flow.  In 
particular,  the  effect  of  the  rotor  blade  wake,  blade  shank  wake,  and  hub  wake 
on  the  separation  of  the  flow  over  a lifting  helicopter  in  forward  flight  was  analyzed. 

Specific  capabilities  of  Program  SHAPES  include; 

• Complete  3-D  representation  of  fuselage/pylon/rotor/hub. 

• Capability  to  input  arbitrary  aerodynamic  paneling  geometry.  i 

• Preprocessor  for  simplified  helicopter  geometry  to  enhance 
parameter  studies. 

• Capability  to  run  all  or  part  of  helicopter  configuration. 

• Capability  to  analyze  separation  effects  on  fuselage  and  pylon. 

• Automatic  generation  of  rotor  and  hub  wake  dynamics. 

• Capability  to  run  multibloded  rotor  problem. 
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The  method  is  very  flexible  and  simple  to  use;  the  use  of  quadrilateral 
hyperbolaidal  elements  (which  can  be  used  to  yield  any  arbitrary  closed  surface),  de- 
fined in  terms  of  their  corner  points,  is  one  of  the  original  features  of  the  method. 
Other  features  of  the  method  are  the  simplicity  of  the  expression  for  the  coefficients. 
This  makes  the  method  extremely  efficient  from  a practical  point  of  view.  Also,  the 
method  is  accurate  and  relatively  fast,  despite  the  fact  that  no  effort  has  been  made 
to  minimize  the  computational  time. 

In  contrast  to  existing  methods,  which  in  many  instances  require  an  exten- 
sive user's  backgraund  in  aerodynamics  and  familiarity  with  the  specific  method,  the 
present  code  requires  very  limited  human  intervention  and  is  extremely  easy  to  use. 

It  is  emphasized  that  the  only  required  inputs  to  the  program  are  the  locations  of  the 
corner  points  of  the  aerodynamic  elements  and  certain  specifications  to  the  problem; 
i.e.,  Mach  number,  rotor  angular  speed,  angle  of  attack,  and  the  location  of  the 
separation  line.  In  addition,  the  geometry  of  the  wake  is  automatically  generated. 
Present  outputs  available  from  the  computer  program  include  the  aerodynamic  coeffi- 
cients, i.e.,  lift  and  induced  drag,  as  well  as  the  pressure  and  velocity  flow  field 
around  the  configuration. 

It  should  be  noted  that  the  CPU  time  for  the  evaluation  of  the  coefficient 
matrix  is  a large  percentage  (typically  80  to  90  percent)  of  the  total  CPU  time.  An 
improved  computational  scheme  (Reference  50)  has  recently  been  developed.  This 
scheme  consists  of  combined  numerical  (Gaussian  quadrature)  and  analytical  pro- 
cedures for  the  evaluation  of  the  source  and  doublet  integrals  used  in  the  program. 

The  method  employed  in  SOUSSA  (Reference  50)  is  very  accurate.  However,  the 
analytical  evaluation  of  the  coefficients  (which  is  necessary  for  "near  elements")  re- 
quires an  unnecessary  amount  of  CPU  time.  A hybrid  numerical-analytical  scheme 
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for  the  evaluation  of  the  coefficient  is  used.  This  scheme  yields  on  appreciable  re- 
duction of  computer  time  (80  to  90  percent  of  the  original  time;  i.e.,  after  this 
modification,  SOUSSA  requires  one-tenth  of  the  original  CPU  time  for  the  evaluation 
of  the  coefficients).  Larger  savings  are  expected  for  complex  configurations,  since 
in  those  cases  a large  number  of  elements  is  required  simply  to  describe  the  geometry 
accurately.  Hence,  it  is  recommended  that  the  optimal  combination  of  the  various 
schemes  as  well  as  the  trade-off  between  accuracy  and  computer  time  be  investigated. 

Comparison  with  the  results  of  References  44  and  46  - 49  indicates  that  the 
present  results  are  in  excellent  agreement  with  existing  ones  and  demonstrate  the  flexi- 
bility and  usefulness  of  the  program.  A number  of  improvements  are  now  under  in- 
vestigation. It  is  seen  that  the  present  wake  model  does  not  take  into  account  the 
concentration  of  the  vortex  sheet  into  a tip  vortex.  This  tip  vortex  concentration 
has  a strong  effect  on  the  flow  field  around  the  helicopter  fuselage/rotor  configuration. 

Also,  a more  complete  separation  representation  can  be  incorporated.  Improvement 
can  be  accomplished  by  moving  the  separation  line  and  by  varying  the  strength  of  the 
isolated  vortex  element  to  verify  the  area  of  separation.  It  is  recommended  that  the 
aforementioned  improvements  be  included  in  Program  SHAPES. 

The  present  method  has  potential  application  in  the  design  of  helicopter 
configurations  because  it  provides  an  analytical  capability  to  the  engineer  that  is 
simple  to  apply  and  which  can  be  used  to  develop  a low-drag  profile  as  well  as  to 
explore  problem  areas. 
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i i APPENDIX  A 

I SOURCE,  DOUBLET,  AND  VORTEX  LAYERS 

’ I 

In  this  appendix  the  main  properties  of  source,  doublet,  and  vortex  layers  ore 
considered.  The  equivalence  of  vortex  layers  and  doublet  layers  is  also  considered. 

A.l  SOURCE  LAYER 

A distribution  of  sources  over  a surface  T is  called  o source  layer.  The  po- 
tential due  to  a source  layer  Is  given  by 

J/  'S4^  '*■" 

i 

where 

r = I P*  - P I (A-2) 

■ j 

! and  I j is  the  intensity  of  the  source  distribution.  Note  that  a source  layer  does  not  yield 

; , any  discontinuity  on  «p  therefore,  the  tangential  components  of  the  velocity  are  con- 

tinuous on  v.  However,  the  normal  component  of  the  velocity  due  to  a source  layer, 

i ^ , is  discontinuous.  In  order  to  determine  the  discontinuity  of  consider  the  nor- 

I n ' 


mal  derivative  of  Equation  (A-1)  at  a point  P^  onT.  This  yields 
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Next  consider  on  the  surface  E a small  circular  surface  element,  T ^ , with  center 
and  radius  e.  Then  isolate  the  contribution  of  Eg  , 


Note  that  when  e is  sufficiently  small,  E^  is  approximately  planar.  Then  considering  a 
frame  of  reference  X,  Y,  Z with  origin  at  P*  and  the  Z-axis  directed  like  the  normal 


where  (see  Figure  A-2) 


9 = tan"'  (Y/X) 


and  Ij  is  the  value  of  at  . 
o 

Note  that 


{A-6) 


(A-7) 


where  the  upper  (lower)  sign  holds  if  approaches  from  Z*  ^ 0 (Z*  < 0).  There- 
fore,  the  discontinuity  in  on  T is  given  by 


?n 


= 'S. 


(A-8) 


wi  th 


3cpj 
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I Sn  . 


dn  i 


(A-9) 


where  1 indicates  the  side  for  positive  Z,  i.e.,  the  side  toward  which  n is  pointing, 
and  2 indicates  the  opposite  side.  For  a horizontal  surface  A I I “1  lu 
_ jf  n j,  the  upper  normal,  where  u and  t indicate  upper  surface  and  lower 

surface,  respectively. 

The  above  result  could  be  anticipated  by  thinking  of  a uniform  source  distri- 
bution on  any  infinite  plane.  The  flux  per  unit  area  emanating  from  the  layer  is  Ij  . 
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TK®  flux  ip®r  unit  or®a^  011  each  side  of  the  surface  is  2;  therefor®,  the  flow  field 
is  a uniform  velocity  of  absolute  value  Ij^  2 away  from  the  loyer.  The  value  of  -^  'n 
is  I5  2 o»^  the  side  of  n and  -l^.  2 on  H\e  other  side.  Therefore,  the  discontinuity 
is  equal  to  l<.  Of  evaluoted  in  the  direction  of  h^. 


A. 2 OOllBLFT  LAYFRS 


A suifvsce  distribution  of  dcHibleN,  \ , \sheie  n i>  the  luvmol  tv'  the 

v>n  V ' 

surface  1'  ^pointitrq  from  sivie  1 h''  sivfe  2^  is  collevi  vi  vk'ublet  Iviyer  . The  pv'tentivil  vfi'e 


to  o doi'blet  loyei  is 


.0' i (ft) 


tA-IOI 


where  is  the  intensity  of  the  doublet  distribution.  Note  thot  vi  vK'ublet  lovei  slelds 
o discsvrtinuity  on  cy'^.  In  cvdet  to  determine  this  disvv'ntinuity , Vv'nsidei  the  wime 
process  usevT  for  o source  element.  |ntrcHfucii\vj  the  surface  vfefined  in  Subsectiv'n 
A,  1 , Equofiott  (A-IOT  may  be  rewritten  as 

II  'll  — ' .1.1  '0  - 

^ v-s-  ^1'  \4"i/  .'n  \4”i  ’ 

t 

Next  note  thvst  vissr>minvi  I,.,  CvvrsNsnt  wittrin  vrnvt  eviiivsl  K'  li'^  , vsne  vsl'Kiins 

■ l»  ‘ V' 

(see  EqiKition  (A-5TT 

II  i ^IzLViT  Ip  11  n -vf— ) dXvtN 

f .'i>V4”r/  cs  \4”i  ' 


Ip  ^ f - - R-'R 

b 4.-r'^ 
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Therefore,  the  discootinuity  i$ 


° ''d 


I 1 / ^ RdR 

° ^ 0 (r2  t Z?  ^ 


I 1 r 
2 k^Tz? 


2 Z. 


(A-12) 


(A-13) 


SOLID  ANGLES 


It  will  be  useful  later  to  interpret  doublet  distributions  in  ternrs  of  solid 


angles.  Note  that  (see  Figure  A-3) 


= c„  X £ 

I /■/'.  din 

'^JJ  D T 

In  particular,  for  a closed  surface  and  Ip  = constant.  Equation  (A-14)  becomes 


(A-14) 
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(A-15) 


f 


According  to  Equation  (A-22),  Equation  (A-23)  may  be  rewritten  as 


(A-25) 


where  T is  any  arbitrary  surface  bounded  by  the  vortex  line  C.  Noting  that  j = 0, 
Equation  (A-25)  may  bo  rewritten  as 


(A-26) 


(A-27) 


Equation  (A-26)  indicates  that  the  velocity  v^  due  to  a single  closed  vortex  may  be 
expressed  in  term  of  the  potential,  cpp  , due  to  a doublet  distribution  of  constant 
intensity  equal  to  the  vortex  intensity  F,  over  a surface  E bounded  by  the  vortex 
line  C but  otherwise  arbitrary. 

A. 6 EQUIVALENCE  OF  DOUBLET  AND  VORTEX  LAYERS 

The  above  analysis  may  be  extended  to  a vortex  layer,  since  a vortex  layer 
may  be  thought  of  as  a superposition  of  single,  closed  vortices.  The  extension  of  the 


t 

f 


analysU  is  briefly  presented  in  the  following.  Note  that  (see  Figure  A-4)  using 
Equation  (A-22),  the  velocity  due  to  a vortex  layer  may  be  written  as 


now  if 

dl  = d^^a^ 

then  Equation  (A-28)  becomes 


(A -28) 


(A -29) 


Introducing 


which  shows  the  equivalence  of  a vortex  layer  to  a doublet  layer  of  intensity  D,  Note 
that 


(A-34) 


which  is  the  relation  necessary  to  determine  the  vorticity  of  the  vortex  layer  from  the 
intensity  of  the  doublet  layer. 
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arbitrary  vectors 

base  vectors  (see  Equation  (121)) 

base  vectors;  0 = 1,3 

= a • a - ; covariant  component  of  the  metric  tensor 
o p 

adjoint  base  vectors;  o = 1 , 3 

3 Y 

controvoriont  component  of  the  metric  tensor  where  ® ~ ^ 

arbitrary  constants 
source  coefficient 

source  integral  (see  Equation  (107)) 

blade  chord 
arbitrary  constants 
drag  coefficient 


doublet  coefficient 


lift  coefficient 


doublet  integral  (see  Equation  (107)) 
pressure  coefficient 

th 

pressure  coefficient  at  centroid  of  h element 

thrust  coefficient  (see  Equation  (8)) 

domain  function  (see  Equation  (55)) 
arbitrary  functions 
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LIST  OF  SYMBOLS  (Continued) 


Ppp'^pm'^mp' 


= 0;  wake  coefficient 

woke  integral  (see  Equation  (107)) 
adjoint  base  vectors;  0 = 1,3 
adjoint  base  vectors;  p = 1,  3 


~ 9p  * Oy*  covariant  component  of  the  metric  tensor 

contra  variant  component  of  the  metric  tensor  where  g*^  ~ 

rotor  incidence  angle  with  respect  to  freestream  velocity 

unit  vectors  in  x,  y,  z directions,  respectively 

length  of  isolated  vortex  element 

(subscript)  denotes  lower  surface 

vortex  line 

number  of  elements  on  the  woke 

direction  normal  to  the  surface  of  the  body 

normal  to  surface  of  the  body  at  point  P 

number  of  elements  on  the  surface  of  fuselage  and  rotor 

number  of  elements  along  the  circumferential  direction  of 
one  woke  spiral 

number  of  wake  spirals 
pressure 

freestream  pressure 


location  of  the  corner  points  of  the  hyperboloidal  element 


t 


LIST  OF  SYMBOLS  (Continued) 


P0»Pi»P2'P3 

P 

P(0 


Ph(0 

pTE 

m 

Pw 

P* 

Q 

s 

Q 


R 

R, 


mn 


parameters  describing  the  geometry  of  hyperboloidal  elements 
(see  Equation  (119)) 

point  on  surface  of  body 

dummy  point  on  surface  Z 

centroid  of  elements  on  surface  of  fuselage  and  rotor 

centroid  of  element  ^^(t) 

centroid  of  m*^^  elements  on  the  wake 


corner  points  along  the  trailing  edge  of  m^^  wake  element 

point  on  the  wake 

control  point  in  the  flow  field 

location  of  separation  line 

location  of  isolated  vortex  line  that  is  used  to  represent  the 
presence  of  vorticity  in  the  separated  wake 

radial  distance  (see  Equation  (43)) 

vector  denoting  the  distance  between  control  point  and 
dummy  point 

rotor  radius 

radius  of  ^2  (*®®  Figure  5) 

function  that  relates  the  value  of  cp  at  the  trailing  edge  with  the 
values  of  ^ at  the  centroid  of  the  element  (see  Equation  (1 12)) 


time 


time  when  point  on  the  wake  is  at  x 


M 
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V 

V 


V ,v  ,v 
x'  y'  z 


V 

V 


V 


V ,V  ,V 
x'  y'  z 
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thrusf 

(subscript)  denotes  upper  surface 

induce  rotor  velocity  (associated  with  the  convection  of  the 
rotor  woke) 

freestreom  velocity 

velocity  induced  by  the  rotor  (see  Equation  (5)) 
perturbation  velocity 

P®*'turbotion  velocity  in  x,  y,  z direction,  respectively 

covariant  components  of  vector  v 

velocity  induced  by  isolated  vortex 

flight  velocity 
velocity  vector 

volume  surrounding  the  surface  of  the  body  and  woke 
(see  Figure  5) 

resultant  flow  velocity  in  the  far  wake 
average  value  of  velocity 

velocity  of  body 

components  of  velocity  in  x,  y,  z directions,  respectively 

Cartesian  coordinates 
location  of  a point  on  the  wake 

vector  describing  the  vortex  lines 

spanwise  coordinate  measured  from  midspan  of  an  elliptical 
monoplane  wing 

blade  root  cut-out 
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LIST  OF  SYMBOLS  (Continued) 


vortex  intensity 

midspan  value  of  elliptical  circulation  distribution  of  a 
monoplane  wing 

n a = Y • 

Kronecker  delta  = j q a / Y 
radius  of  surface  (see  Figure  5) 

= Vx  v (see  Equation  (75)) 

controvarient  components  of  vector  Q}  ot  = 1,3 

rotor  angular  position 

parameter  describing  vortex  lines 

surface  which  surrounds  the  volume  V 

surface  of  the  body 

branch  surface  of  isolated  vortex  line  (see  Figure  8) 
surface  of  the  fuselage 
surface  of  the  fuselage  wake 
h^  element  on  the  surface  of  the  body 

i*^^  elements  on  the  surface  of  the  fuselage,  rotor,  and  wake 
m*^^  elements  on  the  wake 

elements  on  surface  of  fuselage  and  rotor 
surface  of  the  rotor 
surface  of  the  rotor  wake 
surface  of  the  wake 
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E' 

w 


m 

9 

V 

9 

'Ph 

9, 


9, 


m 

TE 

m 


surface  surrounding  each  wake 

surface  surrounding  the  body  and  wakes  (see  Figure  5) 

spherical  surface  of  radius  e and  center  P*  (see  Figure  5) 

spherical  surface  of  radius  R_  and  center  P*  (see  Figure  5) 

2 

time  necessary  for  material  point  of  wake  to  move  from  p|^^  to 
p'^  with  velocity 

perturbation  potential  (see  Equation  (127)) 

velocity  potential  in  absence  of  the  isolated  vortex 
(see  Equation  62) 

th 

value  of  potential  at  centroid  of  h element  on  the  surface  of 
fuselage  and  rotor 

th 

value  of  potential  at  centroid  of  m^  element  on  the  wake 
TE 


value  of  9 at  p 


m 


th 


9. 


value  of  potentiol  at  centroid  of  n elements 


'Ppp'^pm'^Pfiip' 

*Pmm 

9q/9] »92'93 

f 
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value  of  the  potential  at  nodal  points  of  the  element 

parameter  that  yields  a continuous  potential  distribution  on 
the  element 

velocity  potential  (see  Equation  (17)) 

inflow  ratio  (see  Equation  (7)) 

advance  ratio  (see  Equation  (6)) 

surface  coordinate  associated  with  the  hyperboloidal 
element  (see  Figure  10) 


I 


a 


generalized  curvilinear  coordinates;  a = 1,  2,  3 
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P 

X 


Q 

V 


air  density 

wake  deflection  angle  of  the  rotor  wake  geometry  with 
respect  to  the  freestream  (see  Figure  3) 

downwash  = ^ 
dn 

normal  wash  at  centroids  of  elements  S 

n 

angular  velocity  of  rotor 

Laplacian  operator 

gradient  operator  = i—  + 

9x  9y  5z 
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